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Abstract 

Several current projects aim at building a large water-Cherenkov detector, with a fiducial volume 
about 20 times larger than in the current Super-Kamiokande experiment. These projects include 
the Underground nucleon decay and Neutrino Observatory (UNO) in the Henderson Mine (Col- 
orado), the Hyper-Kamiokande (HK) detector in the Tochibora Mine (Japan), and the MEgaton 
class PHYSics (MEMPHYS) detector in the Frejus site (Europe). We study the physics potential 
of a reference next-generation detector (0.4 Mton of fiducial mass) in providing information on su- 
pernova neutrino flavor transitions with unprecedented statistics. After discussing the ingredients 
of our calculations, we compute neutrino event rates from inverse beta decay {i>eP — > e~^n), elastic 
scattering on electrons, and scattering on oxygen, with emphasis on their time spectra, which may 
encode combined information on neutrino oscillation parameters and on supernova forward (and 
possibly reverse) shock waves. In particular, we show that an appropriate ratio of low-to-high en- 
ergy events can faithfully monitor the time evolution of the neutrino crossing probability along the 
shock-wave profile. We also discuss some background issues related to the detection of supernova 
relic neutrinos, with and without the addition of gadolinium. 

PACS numbers: 14.60.Pq, 97.60.Bw, 29.40.Ka 
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I. INTRODUCTION 



The successful operation and the great scientific impact of the Super-Kamiokande (SK) 
water- Cherenkov experiment 0, |^ have motived several research groups to investigate in 
detail the feasibility of a Megaton-class detector |^ for nucleon decay and neutrino physics. 
Three main projects are currently being pursued: the Underground nucleon decay and 
Neutrino Observatory (UNO) project in the Henderson Mine (Colorado, U.S.) 0, the 
Hyper-Kamiokande (HK) project in the Tochibora mine (Japan) |6i], and the MEgaton 
class PHYSics (MEMPHYS) project in the Frejus site (Europe) [7|. These next-generation 
detectors are characterized by a prospective fiducial mass of 0.4 Mton or higher, which will 
allow studies of neutrinos of astrophysical or terrestrial origin with unprecedented statistics 
and sensitivity 0]. 

In particular, a new window would be opened on supernova neutrino physics: by naively 
rescaling the detector mass, it turns out that a 0.4 Mton experiment would observe a 
SN1987A-like signal with a statistics ~ 200 times higher than in Kamiokande (2.14 kton) 
[3], and could promote the current SK upper limit on supernova relic neutrinos P into 
a positive detection. In general, the observation of a large number of supernova neutrino 
events will allow statistically significant spectral analyses in the energy, angular, and time 
domain. In this context, the identification of relatively model- independent spectral features 
is important to disentangle, as far as possible, information related to supernova physics and 
neutrino emission from those related to neutrino fiavor transitions. Recent spectral studies 
which refer to prospective Mton-class water-Cherenkov detectors include: analyses in the 
energy domain to probe neutrino mass- mix ing pa rameters [lol |. to perform multiparameter 
fits including neutrino emission parameters 11, and to identify Earth matter effects 



analyses in the time domain to identify the neutronization burst 0, or signatures of 
shock- wa ve p ropagation effects [10, 16]; analyses in the angular domain to achieve supernova 
pointing jl^; background reduction projects using gadolinium 3, 



In this work, we aim at studying some time dependent and independent observables 
related to (extra) galactic core collapse supernovae, assuming a prospective 0.4 Mton fiducial 
mass water-Cherenkov detector. In particular, we investigate the effects of the supernova 
shock propagation on the observable neutrino signal. We analyze the imprint of the shock 
wave on the time spectra of inverse beta decay events, and show that the ratio between 
the number of events in two suitably chosen energy ranges can actually monitor the time 
dependence of the neutrino crossing probability Ph, thus opening a unique opportunity to 
study shock-induced flavor transitions in "real time". Moreover, we also study the shock- 
wave imprint on the time spectra of events coming from scattering of neutrinos on oxygen 
and on electrons. We point out that the shock wave produces a characteristic distorsion 
of the oxygen event spectra in time domain. The elastic scattering spectrum is instead 
rather insensitive to the shock-wave propagation and might be thus used to track the overall 
decrease of neutrino luminosity. 

We also study other observables that are basically insensitive to shock-wave effects. In 
particular, we discuss the sensitivity of a 0.4 Mton detector in detecting total event rates for 
various interaction processes, supernova relic neutrinos and neutrinos from silicon burning 
before core collapse. 



^ We mention that large liquid argon detectors (not considered in this work) can also provide important 
information on supernova neutrinos; see. e.g., the recent works |2fi l2lL l22l 12,11 124| . 
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The structure of this paper is as follows. In Sec. II we present the ingredients of our 
calculations. We discuss our choice for several input (supernova simulations, neutrino os- 
cillation parameters, effective detection cross sections) needed to compute the SN v signal. 
In Sec. Ill we discuss the results of our calculations for observables which are not sensitive 
to shock waves, i.e., the total number of events in the pre-supernova phase and after super- 
nova explosion, and the energy spectra of the diffuse SN relic neutrino background with and 
without gadolinium loading. In Sec. IV we investigate in detail the shock-wave signature on 
the neutrino time spectra, and present a possible strategy to monitor the time evolution of 
the neutrino crossing probability along the shock wave profile, by means of an appropriate 
ratio of low-to-high energy events. Conclusions and prospects are given in Sec. V. 



II. OVERVIEW OF CALCULATIONS 

In this Section we describe the main aspects and ingredients of our calculations of super- 
nova neutrino event rates. We remind that, in general, numerical simulations of supernova 
explosions provide the unoscillated double differential neutrino distribution in energy and 
time, 

F° = ^ (1) 

dEdt ' ^ ^ 

where z/ = z/g, z7e and in standard notation (25l.l2fil| [x indicating any non-electron flavor). 
Such initial distributions are in general modified by flavor transitions (see, e.g., [l^ 126^). 

Fl-^F^ , (2) 

and must be convoluted with the differential interaction cross section Of. for electron or 
positron production, as well as with the detector resolution function i?e, and the efficiency 
e, in order to finally get observable event rates [lo| . 

Ne = Fy®a^®R^®E . (3) 



A. Neutrino spectra in energy and time 

To our knowledge, only the Lawrence Livermore (LL) group has published in detail [2j 
a successful simulation of supernova explosion and neutrino emission for a time interval long 
enough (~ 14 s) to cover the phenomenon of shock propagation 0, |2^. Since we are 
particularly interested to shock-wave signatures on neutrino flavor transitions, we assume 
the LL simulation (which refers to a supernova progenitor mass of ~ 20 solar masses) as our 
reference input for t > (time after bounce). 

For t < 0, Odrzywolek, Misiaszek and Kutschera (OMK) [s^ have recently estimated 
that the silicon burning phase preceding supernova explosion can release an energy of about 

5.4 X 10^° erg in neutrino-antineutrino pairs (with : v^J^ 5 : 1) for about a couple of 
days, assuming a 20-solar mass progenitor (see, however, [SlJ for different time estimates). 
In the absence of more detailed information, we assume a uniform neutrino flux from silicon 
burning for two days before t = 0, with emission parameters taken from |30| . 

In both cases (t > and t < 0) we factorize the differential distribution of Eq. ([T]) as: 

^.° = ^^(i^.), (4) 
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FIG. 1: Reference neutrino emission parameters as a function of time for different flavors: z^e, i^e, and v^, 
assuming a progenitor mass of ~ 20 solar masses. Left panels: neutrino emission rate dN^/dt and average 
energy (E) for the silicon burning phase, as taken from the OMK calculations ^20]- Right panels: neutrino 
emission rate and average energy for the supernova neutrino burst, as taken from the Lawrence Livermore 
(LL) group simulation j27| . 



for any flavor (z/ = z/e, and u^), where dN^/dt represents the neutrino emission rate 
(number of v per unit time), while if{E,t) is the normalized (/ dE (p = 1) energy spectrum 
parametrized as in [s^ 



where (E) is the average neutrino energy and a is an energy shape parameter. In general, 
both a and (E) can be function of time. For simplicity, we have taken a = 3 for all flavors 
for t > [13, and a ~ 4.3 (from a fit to the spectrum in [s^) for t < 0. For t > 0, the 
average energies for each flavor are given as a function of time in the LL simulation (27| . 
while for t < they are assumed nearly constant and equal to ~ 1.8 MeV for all flavors in 
the OMK evaluation IsT 
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FIG. 2: Normalized neutrino spectra averaged over time. Left panel: energy spectra in the silicon burning 
phase, as derived from ;30| with shape parameter a = 4.28 [S^. Right panel: energy spectra associated to 
the supernova neutrino burst, averaged in the time interval for shock-wave effects {t > 2), as derived from 
[2^ with shape parameter a — 3 ,32J. 



Figure shows the main characteristics of our reference supernova emission modeL The 
upper and lower panels show the neutrino emission rate dNy/dt and the average neutrino 
energy [E) for different flavors. The left panels refer to the silicon burning phase from the 
OMK calculations jSfl,], while the right panels refer to the supernova neutrino burst from 
the LL group simulation j27|. According to the LL simulation, all flavors have comparable 
emission rates (within a factor ~ 2) in the so-called cooling phase {t ^ 0.5 s). In the 
preceding phase the relative v^. emission rate is higher, and shows a distinct neutronization 
peak at 0.04-0.05 s. For later purposes, we observe that the LL simulations predict limited 
variations (< 20%) of the average neutrino energy in the time range (t > 2 s) relevant for 
shock-wave effects on neutrino flavor transitions. 

Figure El shows the normalized neutrino spectra ip for each flavor, after a flux-weighted 
average over time, for both the silicon burning phase (left panel) and after supernova explo- 
sion (right panel). For later purposes, in the right panel we have restricted the integration 
interval to t > 2 s, i.e., to the interval where shock- wave effects are relevant for flavor tran- 
sitions. Note that, in the silicon burning phase, the spectra are equal for all flavors, and are 
peaked at relatively low energy (< 2 MeV) js^. During the supernova explosion, neutrinos 
have an order of magnitude higher average energy, and are peaked at different energy for 
each flavor [2^ . Of course, the detailed features of the reference emission model in Figs. ^ 
and|21niust be taken with a grain of salt, since the distribution of the total energy in flavor 
and time is currently subject to large uncertainties, which may be reduced in more advanced 
future simulations 13311 . 
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B. Neutrino crossing probability and shock waves 



Assuming standard three-neutrino mixing, the parameters relevant to flavor transitions in 
supernovae are the two independent squared mass difference 5m? and Aw? and the mixing 
angles 6^12 and 6'i3 (see, e.g. , 33|). For 5rn? ~ 8 x 10^^ eV^, as currently indicated by reactor 
and solar neutrino data [35[ , the dependence on dm? actually vanishes, while the dependence 
on Am^ and on ^13 can be essentially embedded in the so-called crossing probability Ph (up 
to Earth matter effects that, for simplicity, we do not consider in this work). In general, Ph 
takes the same form for both neutrinos and antineutrinos 13^ 



Ph = PHiAmyE,sm^en,Vix,t)) , (6) 

where V is the neutrino potential profile js^ at radius x and time t, 

Vix,t) = V2GFNeix,t), (7) 

and Ne is the electron density. When needed, we fix Am^ ~ 2.4 x 10~^ eV^ from atmo- 
spheric and accelerator data js^l, sin^6'i2 — 0.3 from reactor and solar data js^, and take 
representative value of sin^ ^13 below current upper bounds [s^- 

Matter effects are potentially relevant when the potential equals the neutrino wavenumber 
kn = Am^/2E, 

V{x,t)c^kH. (8) 

Indeed, it has been shown in fsi'] that, even for non-monotonic supernova density profiles, 
the evaluation of Ph in Eq. (jU)) can be obtained with good approximation through a simple 
procedure involving the ordered product of matrices embedding the local crossing probabili- 
ties ji^l, evaluated at all points where Eq. (jHJ is fulfilled. This approximation is particularly 
useful to study the effect of non trivial density profiles, such as those induced by shock 
waves. 

In a seminal paper j2^, Schirato and Fuller noted that the propagation of the forward 
shock wave in the LL simulation [2^ could influence neutrino flavor transitions a few seconds 
after core bounce. In addition to forward shock effects, it has been recently pointed out 
that a second ("reverse") shock front, is also expected to propagate behind the forward 
one (at a lower velocity), although a detailed description is not yet possible within current 
numerical experiments. 

Figure El shows the simplified shock-wave profiles used in this work. The upper panel 
refers to the neutrino potential V in the presence of forward shock only j2l|, using the 
same parametrization as in [39]. The main features of the forward shock profile are a sharp 
discontinuity at the shock front (which can induce a strongly nonadiabatic transition) leaving 
behind an extended rarefaction zone. In the lower panel we have graphically adapted the 
results of the simulation in to account for a reverse shock, characterized by a smaller 
discontinuity at the front. In both panels, we also show the band spanned by the neutrino 
wavenumber kn = Air? /2E for E G [2,60] MeV. The condition for large matter effects 
{y ~ kn) implies then that at relatively early (late) times only the static profile (the 
rarefaction zone) is relevant, while for intermediate times (e.g., t =2, 4 or 8 s in Fig. EI) the 
propagation of the shock front (s) must be accounted for.^ 



^ Neutrinos from silicon burning may also experience flavor transitions along the static profile (not consid- 
ered in li^)- In fact, for a typical energy £^ ~ 2 MeV (see Fig. |21 left panel), the condition V{x) ~ ku is 
fulfilled at a; ~ 10"* km for the static profile, i.e., below the estimated sihcon core radius (xsi ~ 4 x 10'^ km). 
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FIG. 3: Radial profile of the neutrino potential V(x) at different times t (1, 2, 4, 8, and 16 s). Upper panel: 
Our simplified profile for the case of forward shock (see also Lower panel: Our simplified profile 

for the case of forward plus reverse shock, adapted from In both cases, the static profile (t < s) is 
also shown. The band within dashed lines marks the region where matter effects are potentially important 
{V ~ kn for E = 2-60 MeV). 



Figure ID shows our calculation of Ph as a function of time for four representative values 
of sin^ 6'i3, and for a relatively high neutrino energy {E = 45 MeV, useful for later purposes). 
The left panels refers to the case with forward shock only, where the strong variations induced 
by the passage of the front discontinuity and then by the rarefaction zone are marked by 
vertical dotted lines. We refer the reader to [3^ for a thorough discussion of forward-shock 
effects on Ph- The right panels refer to the case with forward plus reverse shock, where 
the effects of the reverse shock appear to change dramatically the crossing probability at 
intermediate times. Qualitatively, the nonadiabatic transition at the reverse shock front can 
partially "undo" the effect of the analogous transition at the forward shock front; this is 
particularly evident, e.g., in the two upper panels. We refer the reader to ,116] for further 
discussions about the time dependence of Ph in the presence of a reverse shock. Despite 
its complexity, such time dependence might be monitored surprisingly well in future water- 
Cherenkov detectors, as we shall see in Sec. IV. 
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FIG. 4: Crossing probability Ph as a function of time, at relatively high neutrino energy {E — 45 MeV). 
Left panels: case with forward shock only. Right panels: case with forward plus reverse shock. From 
top to bottom, the mixing parameter sin^ 6'i3 takes the values 10~^, 10~^, 10"'*, and 10~^. The function 
Pnit) changes rapidly at the times indicated by dotted vertical lines, i.e., when the static profile is first 
perturbed by the forward shock front (followed by the reverse shock front in the right panels) and then 
by the rarefaction zone. Notice how the reverse shock partially "undoes" the Ph variation induced by the 
forward shock. 



C. Neutrino spectra and "critical energy" 

The neutrino fluxes at the supernova exit {F„) are hnear combinations of the initial fluxes 
with coefficients governed by Ph and by the mixing angle 612, as well as by the mass 
spectrum hierarchy (normal or inverted). In particular, for normal mass hierarchy it is (see, 

e.g., Ei) 

F^^ ~ cos^ ^12^0^ + sin^ ^12^0^ , (9) 
F,^ ~ sin^ ^i2PhF° + (1 - sin^ 9i2Ph)FI , (10) 

while, for inverted hierarchy, 

^ cos^ O^^PhFI + (1 - cos^ e^2PH)Fl , (11) 
F,, ~ sin^ ^i2i^° + cos^ e^^Fl . (12) 
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Notice that, for I/g, the case of inverted hierarchy with Ph = 1 is indistinguishable from 
the case of normal hierarchy; similarly, for Ug, the case of normal hierarchy with Ph = 1 is 
indistinguishable from the case of inverted hierarchy. 

For later purposes, it is useful to introduce the concept of "critical energy" Ec defined 
as the energy where the initial z7g and fluxes are approximately equal, 

(13) 

At the critical energy, flavor transitions are not effective in the antineutrino channel, and 
the flux at supernova exit equals the initial one, 

F^SEc)^FliE,) (14) 

for both normal and inverted hierarchy. 

In the LL simulation [2^, the critical energy is relatively well defined and stable for 
t > 2 s (the range relevant for shock- wave effects); in this time interval, the Ve and 
emission rates are comparable, and their average energies do not vary too much in time (see 
Fig. Q). Within the LL simulation, the critical energy is then approximately defined by the 
crossing point of the time-averaged Vf. and spectra in Fig. |2l i.e., ~ 20 MeV. This 
energy provides a sort of "no oscillation benchmark", in comparison with higher energies 
(say, E ~ 40-50 MeV), where ^ F-^^ (see Fig. ^ and fiavor oscillation effects do not 
cancel. Of course, the next real supernova explosion might be characterized by a "critical 
energy" different from this or other simulations.'^ We shall discuss in Sec. IV A a possible 
way to circumvent our "a priori" ignorance of the real value of Ec, when seeking signatures 
of shock-wave effects. 



D. Neutrino interactions and detection 

In large water- Cherenkov detectors, the interaction processes which can provide observ- 
able supernova neutrino event rates are: (i) inverse beta decay (possibly with neutron capture 
in Gd tlSj); (ii) (anti) neutrino scattering on oxygen; and (iii) (anti)neutrino scattering on 
electrons. For the process (i), we take the differential cross section from ^l| (we recall that 
the positron energy Epos closely tracks the neutrino energy E |4l[). For the process (ii), we 
take the total cross section from [i^], where we have assumed the same energy distribution 
for all the channels. The differential cross sections for the processes in (iii) are well-known 
(see, e.g., |T3). 

For any interaction process, we fold the differential cross sections for production with 
a Gaussian energy resolution function of width A, and apply a sharp cut to events with 
measured energy below a threshold value i^thresh- The value of A is determined by the 
photocatode coverage. For a coverage comparable to the current one in Super-Kamiokande II 
(~ 20% of the area |43|), as envisaged in typical Mton-class projects |8i], we derive from the 
information in [i^ that 

A/MeV ~ 0.6^E/MeV . (15) 



^ A remark is in order. The concept of "critical energy" is well defined only if the Pe and Vx energy spectra 
turn out to have different shapes and thus a crossing energy point, as is the case for the LL simulation l27j 
at Ec — 20 MeV. Recent simulations [s^ predict 9^ and Vx spectrum differences which, although smaller 
than in the LL case, in general still allow a meaningful definition of the critical energy. 
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FIG. 5: Effective neutrino interaction cross sections as a function of energy, including energy resolution 
and threshold effects. See the text for details. 

We also adopt the conservative value -Ethresh = 7 MeV from 0. The detection efficiency 
e is assumed to be 1 above threshold. Unless otherwise noticed, the fiducial volume for 
supernova neutrino detection is assumed to be 0.4 Mton j3,0|- 

Figure El shows the effective cross sections as a function of neutrino energy for the various 
processes discussed before, including resolution and threshold effects. Notice that, by adding 
Gd [l^, the neutron capture signature allows to detect sub-threshold inverse beta decay 
events. 

III. OBSERVABLES NOT PROBING SHOCK WAVES 

In this Section we discuss the results of our calculations for time-integrated observables, 
such as the total number of events from a single supernova explosion or from the diffuse 
supernova neutrino background, which are basically insensitive to shock-wave effects. Oscil- 
lation effects are simply accounted for by varying Ph in its full range [0, 1] for both normal 
and inverted hierarchy (for a static profile^this is equivalent to vary sin^ 6'i3 in the currently 
allowed range, sin^ 6*13 ^ few%, see e.g., |lO|). 

A. Total number of events 

Figure ini shows the total number of events expected in a 0.4 Mton detector for various 
interaction processes, as a function of the supernova distance in kpc. (For the sake of com- 
parison, the vertical axis on the right refers to the number of events for a SK-like volume of 
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FIG. 6: Number of events expected in 0.4 Mton detector (left y-axis) or in a SK-like detector (right y-axis) 
as a function of the supernova distance, for various interaction channels. See the text for details. 
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22.5 kton). Relatively close stars which might evolve into core-collapse supernovae at un- 
predictable future times include Betelgeuse, Mira Ceti, and Antares. For close galaxies 
the rate of core-collapse supernovae may be inferred from their luminosity and morphology 
4^1; representative estimates based on 0, 0| (to be taken within a factor of two) are 



reported in the upper part of the figure for: the galactic center, the Large Magellanic Cloud 
(LMC), and the M31+M33 spiral galaxies. For each interaction process, the uncertainty due 
to our ignorance of the hierarchy and of Ph G [0, 1] is accounted for by the vertical spread 
of the bands. The relatively different spreads reflect different sensitivities to neutrino flavor 
transitions. Let us now discuss the various classes of events from bottom to top. 

Inverse beta decay events from silicon burning have very low (sub-threshold) positron 
energies, and can only be detected through neutron capture by adding gadolinium [l8|, 
provided that they can be statistically distinguished from background fluctuations. At a 
distance of 1 kpc, we estimate ~ 200 to ~ 800 Si burning events over about two days. This 
large variation is due to the large ratio F^^ : F^^ ^5:1 (see Sec lA), which implies Fy^ ~ 0.2- 
0.8 X F^^, where the lowest value is obtained for inverted hierarchy and Ph = 0, while the 
highest value is obtained for either inverted hierarchy with P^ = 1 or for normal hierarchy. 
The estimated signal rate of ~ 100-400 events/day at 1 kpc should be compared with 
the background rate from spallation neutrinos (e.g., 1780 events/day at the Kamioka site), 
reactor neutrinos (e.g., ~ 500 events/day at the Kamioka site), and from minor contributions 
(natural radioactivity). We have taken such backgrounds from ^] and rescaled them to 
a 0.4 Mton mass. Assuming a typical total background rate of ~ 2500 ± 50 events per 
day (statistical errors only), the silicon burning signal should then be seen with a statistical 
significance of 2-8 standard deviations at a reference distance of 1 kpc. Unfortunately, at the 
galactic center (~ 10 kpc) the estimated silicon burning signal would be 100 times smaller 
and thus unobservable. 

There are better prospects to observe the neutronization burst from a galactic supernova 
by means of elastic scattering on electrons, including contributions from all flavors. The 
contribution of antineutrinos can be neglected, being suppressed both by the original flux 
and by the cross section (see Figs.^andEl). The contribution of non-electron neutrino flavors 
can be derived by using Eq. (fTUI) or Eq. (fT^ plus the unitarity condition Fy^ + F^^ + F^^ = 
F^^ + 2F°^ . In general, the dominant contribution comes from Fy^ (with F^';^^^ ~ sin^ 
either in normal hierarchy with Ph = 1 or in inverted hierarchy), except for the case of 
normal hierarchy with Ph — 0, where F^,^ + F,^^ ~ Fj^^ dominates. Therefore, we expect a 
signal variation by a factor ~ sin^ 9 uO'u J o'u^ ~ 2 (where a is the elastic cross section), due 
to our ignorance of the hierarchy and of Ph- Indeed, at 10 kpc we estimate from ~ 12 to 
~ 22 events during the time interval of the neutronization peak (t G [42, 47] ms in the SN 
simulation, see Fig.^). The smallness of this time interval and the forward-peak signature of 
elastic events make this event sample basically background-free. In conclusion, a 0.4 Mton 
detector might observe the neutronization signal from a galactic supernova with a typical 
statistical significance of 3.5-4.7 standard deviations, assuming the emission model in Fig. ^ 
At the distance of the Large Magellanic Cloud, however, the sensitivity drops dramatically 
[0(1) event, see Fig. IH]. 

Let us now consider the (total) contributions from elastic scattering, absorption on oxy- 
gen, and inverse beta decay, integrated over the time interval t e [0, 14] s. For elastic 
scattering, the slanted band in Fig. |B1 is characterized by the smallest vertical spread, since 
the contributions from all flavors tend to partly cancel flavor transition effects. In a sense, 
the elastic scattering event sample (~ 4 x 10^ events at 10 kpc, easily separable through 
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forward-direction cuts) could be used as a reference value to estimate the original, unoscil- 
lated neutrino flux. The cross section for z/g and Vf, absorption on oxygen becomes rapidly 
larger than the elastic scattering cross sections above ~ 20 MeV (see Fig. EJ, where the 
adopted energy spectra are still sizeable (see Fig. |2)). At = 10 kpc, we obtain about 
1.1 X 10'^ events (±20% uncertainty from ignorance of the hierarchy and of Ph)-, with a 
significant contribution from v^. If the case of normal hierarchy, the contribution might 
provide a handle to flavor transitions not observable in the "canonical" Vf, channel (inverse 
beta decay), as we shall discuss in the context of shock- wave signatures. We remind that 
events from (anti) neutrino absorption on oxygen are slightly backward peaked 4^ . 



Finally, the inverse beta decay channel will provide, in general, such a large statistics 
in a 0.4 Mton Cherenkov detector (~ 2 x 10^ events at 10 kpc) that a handful of events 
might be seen even at a distance as large as ~ 1 Mpc. The high statistics available for a 
galactic supernova explosion will allow many possible spectral analysis, examples of which 
will be given in Sec. IV in the context of shock-wave effects. Here we simply notice that, 
since the angular distribution of e"*" from inverse beta decay is slightly forward peaked (see, 
e.g., 01), a statistical separation from oxygen absorption events appears quite feasible, not 
only for the total event sample, but also for subsamples (e.g., time or energy bins). For 
instance, we have checked that, for an O(IO^) inverse beta decay event sample, the slope of 
the positron distribution in cost? (where °d is the scattering angle between z/g and e^) changes 
significantly — with respect to the very small statistical errors — when O(IO^) oxygen events 
are added (not shown); further studies are needed, however, to make quantitative statements 
in this sense. 



B. Supernova relic neutrinos 

In this section we study the flux of neutrinos coming from all past core-collapse supernovae 
(so-called Supernova Relic Neutrinos, SRN). Since the Super-Kamiokande SRN upper limit 
is only a factor of a few above typical flux estimates , the SRN detection through inverse 
beta decay appears a feasible goal in Mton-class detectors, especially if the background 
can be suppressed by capturing the neutron with gadolinium [13|. SRN detection might 
be even closer ji^ in view of recent experimental data from the Galaxy Evolution Explorer 
(GALEX) mission [i^l, which suggest an increase of the estimated star formation rate (SFR) 
and thus of the SRN flux [i^. In the following, we present quantitative estimates of the SRN 
signal (for a 0.4 Mton detector), based on the emission model discussed in Sec. II A and on 
the same astrophysical input as we used in |4^, except that we adopt here the lowest SFR 
function allowed by the (dust-corrected) GALEX and other data (i.e., the curve labeled as 
"min Apuv" in Fig- 5 of [i^), consistently with the "concordance SFR" advocated in ||46|. 
Our background estimates refer to the Kamioka site, for definiteness. 

Figure d shows, in the upper left panel, the absolute relic neutrino spectrum expected in 
the two extreme cases of normal hierarchy and of inverted hierarchy with Pu = 0, the latter 
being the most favorable in term of event rates. The abscissa is labelled by the measured 
positron energy. The vertical scale gives the number of events in 5 MeV bin of the measured 
positron energy, for an exposure of (0.4 Mton) x (4 years); for the sake of clarify, however, are 
shown as continuous curves rather than as 5-MeV bin histograms. In the upper left panel, 
we also show the background from invisible muons (with normalization taken from the SK 
measurement j^, rescaled to 0.4 Mton), and the background from low-energy atmospheric 

and Fg, which we have calculated using the recent FLUKA fluxes ji^. We remind that 
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FIG. 7: Supernova relic neutrino signal and background for a hypothetical 0.4 Mton detector at the 
Kamioka site, without and with gadolinium (left and right panels, respectively). Upper panels: Energy 
spectra of the signal and of atmospheric and reactor neutrino backgrounds. Lower panel: subtracted signal 
with total statistical errors (for an exposure of 4 years). 



in the SRN energy range, the z/g's interact dominantly through inverse beta decay and 
sub dominant ly through absorption on oxygen, while v^-s, dominantly scatter on oxygen. For 
sirnplicity, we have assumed no atmospheric v^, or z7e oscillations; however, as also remarked 
in |48| , we stress that oscillation effects induced by the "solar" squared mass difference bw? 
cannot be totally neglected in this context, since dm^R^/E ~ 0(1) in the energy range of 
Fig. [7| The reactor neutrino background has been estimated on the basis of the spectrum 
shape in l50| (appropriate to parametrize the high-energy tail of the reactor spectrum), 
with total normalization provided by the observed flux in the KamLAND experiment ^35j |. 
Finally, a ~ 20 MeV hard cut to reject spallation events is also shown. In the upper right 
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panel (as compared to the upper left one) we assume that inverse beta decay events can 
be tagged by neutron capture in gadolinium. In this case, although the reactor background 
cannot be reduced (having the same signature as SRN events), spallation events are fully 
rejected. Concerning atmospheric neutrinos, it is estimated that invisible muon events and 
z/g oxygen absorption events can be reduced by a factor of ~ 5 with gadolinium |0 . 

From the comparison of the upper panels in Fig. [7|we learn that, even with gadolinium, 
the extraction of the SRN signal requires a careful subtraction of (supposedly known) back- 
grounds from the total signal. Therefore, it is important to improve as much as possible 
the knowledge of both the shape and the normalization of the background. In all cases, the 
shape is partly controlled by the energy resolution, which must be known very well; in fact, 
any increase or uncertainty in the resolution width A (due, e.g., to variation in the photo- 
tube coverage) will make some background event "leak" in the SRN sample, especially from 
the reactor neutrino sample, characterized by a steeply falling spectrum. Further studies of 
atmospheric neutrino fluxes 4^ and of their interactions in oxygen i^l in the low-energy 
regime will also be beneficial to the reduce systematic uncertainties in the extraction of the 
SRN signal. 

If we neglect systematics, the SRN number of event Ns (derived from Ntot = Ns + Nb, 
where Nb is the estimated number of background events) is affected by an error ii/A^tot- 
Figure m shows, in the lower panel, the prospective values of Ns±y/N^, for the less favorable 
case of normal hierarchy. The impact of gadolinium is evident both as a reduction of the 
error bars and as a reduction of the energy threshold for SRN detectability. However, even 
without gadolinium, a ~ 2-3a SRN signal could emerge after an exposure of 4 years with 
0.4 Mton. The situation would be more favorable in the case of inverted hierarchy with 
Ph = (not shown in the lower panels of Fig. |2I). 

In conclusion, there are good prospects to reveal a SRN signal with an exposure of a few 
years in a 0.4 Mton detector, especially (but not necessarily) with the addition of gadolinium. 
Further studies are needed, however, for a better characterization of the background. In 
particular, in the SRN energy range, the tails of the reactor and atmospheric spectra, and 
of the energy resolution function, need to be under control to avoid migration of events. 



IV. OBSERVABLES PROBING SHOCK WAVES 

In this Section we study possible signatures of shock-wave effects that might be seen in 
a 0.4 Mton detector. We assume a galactic supernova {d = 10 kpc) and consider first the 
shock signatures in the absolute spectra of inverse beta decay events. We show then that a 
specific spectral ratio can actually monitor the time dependence of the crossing probability 
Ph, thus providing real-time information about the density profile. Finally, we show how 
events from elastic scattering and absorption on oxygen can further help the discrimination 
of shock-wave effects. 



A. Absolute time spectra from inverse beta decay 

Inverse beta decay events can be sensitive to shock-wave effects on only in the case 
of inverted hierarchy [see Eqs. and (jllj) ]. However, inverted hierarchy is a necessary but 
not sufficient condition. For z7g's with E ~ E^, where Ec is the critical energy introduced in 
Sec. II C {Ec — 20 MeV for the our adopted emission model), flavor transition effects will 
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be largely cancelled, including matter effects along the shock profile. No such cancellation 
is expected, however, at higher energy (say, 45 MeV). Therefore, it makes sense to compare 
the expected time spectra of events generated by neutrinos of relatively "low" energy {E ~ 
Ec ~ 20 MeV) and relatively "high" energy {E ~ Eh ~ 45 MeV) for both normal and 
inverted hierarchy, with and without shock in the latter case. 

Figures |H1 and El show absolute time spectra of events for the cases of forward shock only 
and of forward plus reverse shock, respectively, for the same representative values of sin^ ^13 
used in Fig.|3] In both Fig.|Hland|ni the left (right) panels represent the time evolution of the 
number of events within the positron energy bin Epos = 20±5 MeV, (-Epos = 45±5 MeV) for 
the cases of normal hierarchy, inverted hierarchy with static profile, and inverted hierarchy 
with dynamical shock profile. The number of events in the bin Epos = 45 ± 5 MeV is 
generally smaller than for 20 ± 5 MeV. The time spectra in the bin i^pos = 45 ± 5 MeV show 
strong signatures of shock-wave effects for inverted hierarchy, in the form of non-monotonic 
time variations of the event rate, especially at relatively high values of sin^ Oi^ and in the 
case of forward plus reverse shock. The amplitude of such variations is approximately 
contained within the two extreme cases of inverted hierarchy with no shock and of normal 
hierarchy. Conversely, the shock-induced time variations in the bin i^pos = 20 ± 5 MeV are 
significantly smaller, and could be hardly characterized a priori as non-monotonic. This 
feature can provide a sort of "postdiction" of the critical energy Ec, should real shock-wave 
effects be seen in a future supernova explosion: Ec could be identified as the energy where 
nonmonotonic time variations of the observed signal are most suppressed (as compared 
with variations at higher energies). This empirical definition of Ec is sufficient for our 
purposes and, in any case, cannot be made more precise, both because Ec may fluctuate 
with time and because the detector energy resolution function has a nonnegligible width. 
Finally, we observe that the high number of events in each time bin of Figs. |H1 and IHl excludes 
that the spectral variations (or their absence) can be obscured by statistical fluctuations. 



B. Low-to-high energy ratio of time spectra 

The results discussed in the previous section suggest that the signatures of shock wave 
effects can be enhanced by comparing time spectra at the critical energy Ec, where, one has 
by construction 

F^^ {Ec, t) ~ Fl {Ec, t) ~ F° {Ec, t) (16) 

at any t [see Eqs. ()13|) and ()14|) ]. with time spectra at significantly higher energy E^, where 
one expects F^^ -C F^^ (see Figs. [T]and|21) and thus 

Ft,XEh, t) ~ sin2 e^2Fl{EH, t) (17) 

for normal hierarchy [Eq. Q] and 

F^XEu^t) - [1- cos' e,2PH{E,t)]F^XEH,t) (18) 

for inverted hierarchy [Eq. (fTUI) ]. It is convenient to choose Eh so that the flux reduction 
from Ec to Eh [F^^{Eh) / F^^{Ec) <^ 1] is roughly compensated by the increase of the inverse 
beta decay cross section [a{EH) ^ a{Ec)], 

(19) 

Fl{EH)a{EH) • ^ ^ 
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Time spectra of (p, n) e^ (fwd shock) 

Epos = 20±5MeV Epos = 45±5MeV 
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FIG. 8: Absolute time spectra of events from inverse beta decay in a 0.4 Mton detector, in the presence of 
forward shock only, for four representative values of sin^ ^13. The solid, dashed, and dotted histograms refer 
to calculations in inverted hierarchy with shock effects, inverted hierarchy with a static density profile, and 
normal hierarchy, respectively. In each panel, the left (right) panels refer to the energy bin E'pos = 20 ± 5 
(45 ± 5) MeV. 
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Time spectra of (p^ n) e^ (fwd + rev shock) 

Epos = 20±5MeV Epos = 45±5MeV 
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FIG. 9: As in Fig. |51 but for the case of forward plus reverse shock. 
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With such choice for Eh, one can derive from the previous equations the following approxi- 
mate ratio of inverse beta decay events at Ec and Eh'- 



N{E,,t) 
N{EH,t) 



sin^ 9 



1 



12 



(20) 



for normal hierarchy and 



N{E,,t) 1 



(21) 



N{EH,t) l-COS^9uPH{EH,t) 



for inverted hierarchy. 

For our adopted emission model, the condition in Eq. ()19p is fulfilled at Eh ~ 45 MeV, 
which explains a posteriori our choice for the "high" energy bin in Figs. |H1 and IHI Other 
choices for Eh around 45 MeV would simply provide an overall factor on the right hand side 
of Eqs. ()2()|1 and (j2H), with no substantial change in the following discussion. 

Figure El shows the ratio of events in the 20 ± 5 MeV bin with respect to events in 
the 45 ± 5 MeV bin, as compared with the function [1 — cos"^ OuPr {Eh, t)]~^ . Notice that 
this function has the same qualitative behavior as PH{EH,t) itself (see Fig. E}. It can be 
seen that, in inverted hierarchy, shock effects are faithfully monitored by the chosen ratio 
of events, i.e., there is a striking correspondence (even in absolute values, within a factor of 
2 or better) between the solid histograms on the right and the solid curves on the left in 
Fig. up to smearing effects due to the detector resolution. For inverted hierarchy with no 
shock (dashed curves and histograms) there is also a reasonable correspondence of the event 
ratio on the right with the constant value on the left. For normal hierarchy, as expected, any 
information on Ph is lost. Therefore, in the case of inverted hierarchy, the ratio of events 
at "high energy" and "critical energy" appears as a useful tool to track the main variations 
(and possibly the absolute value) of the crossing probability Ph, from which one could get 
precious information about the density gradient along the shock profile. 

Figure [TT] is analogous to Fig. [TUl but for the the case of forward plus reverse shock. 
Despite the more complex structure of the crossing probability function, its behavior is 
faithfully tracked by the ratio of events shown on the right (for inverted hierarchy), except 
perhaps at the lowest value of sin^ 6*13 shown. In any case, the main spikes and valleys of the 
function [1 — cos^ 9 12PH {Eh, t)]~^ are clearly reproduced by the (20±5)/(45±5) MeV event 
ratio, and cannot be mimicked by the very small statistical fluctuations (not shown). Of 
course, relations in Eqs. (fT^. ((201), (PT|) are only approximated, and this is the reason 
why, for example, in the case of normal hierarchy (dotted line in Figs. 1101 and ITT|) or inverted 
hierarchy with static profile and sin^ 6*13 = 10~^ (dashed line in the lower right panels of 
Figs. ^1 and [TT|) the ratio does not appear constant in time. 

Of course, when the next observable supernova will explode, we will not know a priori 
the values of Ec and of Eh fulfilling Eqs. (fTBj) and (fT^ . However, if an imprint of shock 
waves is seen in the absolute spectra through nonmonotonic time variations, E^ can be 
found by seeking the (low) energy range where such variations are suppressed, as already 
mentioned. Concerning Eh, its ignorance is not crucial to track the qualitative behavior of 
PH{t)'- histograms very similar to those shown in the right panels of Figs. ITUl and ITTl would 
be obtained by using, e.g., the 40 ± 5 MeV or 50 ± 5 MeV bin instead of 45 ± 5 MeV, up to 
an overall rescaling of the y-axis (not shown). 
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FIG. 10: Time dependence of the ratio of events between the energy bins -Epos = 20 ± 5 MeV and 
Epos = 45 ± 5 MeV (right), compared with the function [1 - cos'^ 612 Ph {Eh ,t)]-^ with Eh = 45 MeV 
(right). It appears that, for inverted hierarchy, the event ratio on the right tracks rather well the function 
structures on the left (up to smearing effects), thus providing a "shock wave monitor" in real time. 
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FIG. 11: As in Fig. ^| but for the case of forward+reverse shock. 
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C. Time spectra from interactions on oxygen and on electrons 



Both Ue and Vg contribute to oxygen absorption events. Therefore, shock wave imprints 
are expected for both inverted and normal hierarchy. Figure ^1 shows the time spectra of 
such events for both normal and inverted hierarchy, using either the static or the forward 
shock-wave profile for sin^^is = 10^^. It can be seen that, in the presence of the shock, 
the event decrease in time is first steepened and then flattened for both hierarchies; in the 
case of inverted hierarchy, there is even a small dip at 9-10 seconds. Therefore, if the time 
decrease of the event rate is understood from simulations or from independent experimental 
information, the identification of a steepening+flattening behavior in oxygen events might 
be taken as a signal of shock passage (in the absence of alternative explanations). 
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FIG. 13: Time spectra of total elastic scattering events in normal and inverted hierarchy (with forward 
shock profile), for sin^ ^13 = 10~^ and 0. 



Independent experimental information on the time decrease of the (anti)neutrino flux 
can be gained through elastic scattering events. Since neutrinos of all flavors contribute 
to this sample, partial cancellation of oscillation effects occurs, as remarked in Sec. Ill A 
for the total rate in this class of events. Figure El shows that the partial cancellation 
persists also in the time spectra, i.e., there are no large variations between normal and 
inverted hierarchy and with or without crossing probability effects (sin^ ^13 = 10^^ or 0). 
Therefore, the elastic scattering sample is expected to track the overall decrease of flux in 
time, rather independently of flavor transition effects. This "standard candle" will be useful 
for comparison with the oxygen and inverse beta decay sample, where shock-wave effects 
can instead range from moderate to strong. 
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V. SUMMARY AND CONCLUSIONS 



In this work we have analyzed the discovery potential of a 0.4 Mton water-Cherenkov 
detector in providing information on several time dependent and independent observables 
related to (extra) galactic core-collapse supernovae. We have considered elastic scattering on 
electrons as well as inelastic scattering on protons and on oxygen in the detector. 

Concerning time independent observables, we have calculated the total number of events 
as a function of the distance of the exploding star for each detection channel during: a) 
the last pre-supernova (silicon burning) phase; b) the neutronization burst; and c) the core 
collapse phase. We have shown that the low energy {E ^ 5 MeV) thermal neutrinos emitted 
during the Si burning phase a couple of days before the explosion can be identified at 2a 
level over typical backgrounds if the supernova is not too far {d ^ 1 kpc) and the detector is 
loaded with gadolinium. The neutronization burst z/g's can be identified with a significance 
^ 3 standard deviations if (i ^ 10 kpc. From core collapse supernovae located near to the 
galactic center we expect up to O(IO^) events from inverse beta decay events, and O(IO^) 
events coming from elastic scattering and and oxygen absorption events. At the larger 
distance of, e.g., the Andromeda galaxy (M31), we expect 0(10) events, thus opening a 
possible window for near extragalactic supernova surveys. 

We have shown that with the addition of gadolinium, the energy range [10, 30] MeV can 
be fully exploited to detect the supernova relic neutrino flux at the level of several standard 
deviations in a few years. However, even without gadolinium, the SRN background could 
emerge at ~ 2-3a level after an exposure of ~ 4 years, if the atmospheric background is 
under control. 

We have then studied the time spectra of events in the range t G [2, 14] s, in order to 
extract information on the neutrino crossing probability induced by the shock wave propa- 
gation in the stellar envelope during the cooling phase. In particular, we have considered the 
high statistic sample of inverse beta decay events. In this context, it is useful to introduce 
the concept of "critical energy" {Ec — 20 MeV for our reference model), where oscillation 
effect cancel to a large extent. In this way, the time dependence of the rate of the events 
in a bin around the critical energy is mainly affected by the overall decrease of the neu- 
trino luminosity, and is almost independent from the time dependence of the probability. 
We have shown that the ratio between the number of events in a bin close to the critical 
energy {Ef. = 20 ± 5 MeV) and those in a bin centered at a significantly higher energy (e.g.. 
Eh = 45 ± 5 MeV) can faithfully track the neutrino crossing probability (if the mass hier- 
archy is inverted), thus providing a "real-time" movie of shock wave effects into the stellar 
envelope. In particular, the case of forward shock only [28] and of forward-|-reverse shock 
leave significantly different signatures in the time domain, and can thus be discriminated. 

We have completed our work by considering the time dependence of the neutrino event 
rates from oxygen absorption and from elastic scattering. The latter are almost independent 
from shock and flavor transition effects, and can thus be used to monitor the overall decay 
of neutrino luminosity. The former (oxygen) event rate is instead sensitive to the shock 
wave also in the case of normal hierarchy (due to the z/g contribution). If a "steepening- 
fiattening" behavior is observed in the oxygen event time spectrum, it could be interpreted 
as an imprint of the shock wave passage, in both cases of normal and inverted hierarchy. 

In conclusion. Megaton-class detectors offer unprecedented opportunities to study super- 
nova and neutrino properties through high-statistics studies of energy and time spectra of 
supernova neutrino events from different interaction channels. Directions for further studies 
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may include detector simulations of the statistical separation between such channels, more 
refined calculations of the oxygen absorption differential cross section, improved numerical 
simulations of supernova explosion and shock-wave behavior, a better understanding of the 
unoscillated neutrino spectra in time and energy (and of their associated uncertainties), and 
an accurate characterization of the atmospheric and reactor background in supernova relic 
neutrino searches. 
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